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Abstract 

The Higgs boson with a mass Mh ~ 126 GeV has been observed by the ATLAS 
and CMS experiments at the LHC and a total significance of about five standard 
deviations has been reported by both collaborations when the channels — )■ 77 
and H — )• ZZ — t- 41 are combined. Nevertheless, while the rates in the later search 
channel appear to be in accord with those predicted in the Standard Model, there 
seems to be an excess of data in the case of the -ff — )• 77 discovery channel. Before 
invoking new physics contributions to explain this excess in the di-photon Higgs rate, 
one should verify that standard QCD effects cannot account for it. We describe how 
the theoretical uncertainties in the Higgs boson cross section for the main production 
process at the LHC, gg — )• H, which are known to be large, should be incorporated 
in practice. We further show that the discrepancy between the theoretical prediction 
and the measured value of the g'g' — )• — )■ 77 rate, reduces to about one standard 
deviation when the QCD uncertainties are taken into account. 



The Higgs particle [T] has been, at last, observed by the ATLAS and CMS experiments 
at the LHC as a signal with about five standard deviations has been reported by each 
collaboration when the main search channels are combined [2]. This discovery represents 
a triumph for the Standard Model (SM) of particle physics and crowns more than four 
decades of theoretical and experimental endeavour. Now that the Higgs discovery chapter 
is closing, a new and even more challenging chapter is opening: the verification of the 
fundamental properties of the particle and the precise determination of its couplings. This 
program can be started at the LHC since, for the reported mass Mh ~ 126 GeV [2], one 
can have access to the Higgs boson in many production and decay channels [3]. 

At the LHC, the main Higgs production channel is the top and bottom quark loop 
mediated gluon-gluon fusion mechanism gg H: at center of mass energies of ^/s = 7 and 
8 TeV and for a Higgs boson mass of Mh ~ 126 GeV, the inclusive cross sections are about 
(^{gg — -ff) ~ 15 pb and 20 pb, respectively [HEj. The vector boson fusion qq — >■ Hqq 
and the Higgs-strahlung qq HW + HZ mechanisms add only little to these rates, 
respectively ~ 8% and ~ 5.5%, before kinematical cuts are applied [HIS]. The gg ^ H 
cross section is known up to next-to-next-to-leading order (NNLO) in perturbative QCD: 
the ii'-factor defined as the ratio of the higher order to the leading order (LO) |6] cross 
sections is ~ 1.8 at NLO [TJIH] and ~ 2.5 at NNLO [9]. The cross section receives also small 
contributions from the resummation of soft gluons [10] and electroweak corrections [UllIS] . 
Some small corrections that go beyond NNLO accuracy are also available [13] but have 
not been included in the predictions used by the LHC experimental collaborations. It is 
clear that it is this exceptionally large i^-factor that allows for a sensitivity to the Higgs 
boson at the LHC with the presently collected data. The main Higgs search channels 
take advantage of the clean H — )■ 77, H — )• ZZ — )■ ii"^ and H — >■ WW — )■ iuiu final states 
(with i = e,fi), while the significance of other modes such as H^t^t" and VH ^Vbb is 
presently low. The Higgs decay branching ratios are rather precisely known 

The results presented by the ATLAS and CMS Collaborations on their Higgs search [2] 
turn out to be a little surprising. Indeed, while the rates for the H — ZZ search channel 
seem to reasonably agree with the SM expectations, a discrepancy mostly driven by the 
ATLAS results is observed in the channel — )■ 77 which has the largest signal significance. 
Defining the ratios i?xx = '^h^xx/'^h^xx measured cross section in a given search 

channel compared to the theoretical expectation, one finds for the two experiment^ 

ATLAS: i?^^ = 1.90 ± 0.50 , = 1-3 ± 0.6 , 

CMS: i?^^ = 1.56 ± 0.43 , = 0.7 ± 0.5 , (1) 

ATLASeCMS: i?^^ = 1.71 ± 0.33 , = 0.95 ± 0.40 . 

It then seems that the rate in the if — 77 channel is more than two standard deviations 
larger compared to the SM prediction, when the ATLAS and CMS measurements are 
combined. This is a rather exciting situation as, not only the Higgs boson has been finally 
discovered, but in addition it appears to come with hints of some new physics. This could 
be the first (long awaited) signal for beyond the SM physics at the LHC. It is expected 



^The numbers given by the collaborations correspond to the optimal (i.e. which maximises the likelihood 
of the test statistics) value jj, of the signal strength modifier that multiplies the expected cross section such 
that fl = a/a^^. As such, strictly speaking, they are not the true cross section ratios Rxx- We will 
nevertheless assume that they are the same for simplicity and for purposes of illustrating our point. 
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that a large number of studies (including some by the present authors and in addition to 
the many which were done after the first hint for this excess was reported in the 2011 data) 
will be devoted to the explanation of this feature in terms of new phenomenal!. 

However, before doing so, it would be probably wiser to consider more conventional 
explanations for this intriguing excess. The first one would be simply that this is the result 
of a statistical fluctuation (in both signal and backgrounds); after all, many > 2a excesses 
appeared in the recent years and faded away when more data was collectecl!. Another 
possibility would be that the systematical uncertainties in the extremely difficult if — )■ 77 
channel have been underestimatecill; potential experimental problems (if any) could also 
be fixed when more data is accumulated and the detector response better understood. 

A third conventional possibility to explain the i7 — t- 77 excess would be that the 
QCD uncertainties may have been underestimatecH by the experimental collaborations. 
This is the option that we will investigate in the present paper. We will show that if 
the theoretical uncertainties in the prediction of the cross section for the by far dominant 
gg ^ H production process at the LHC are properly included, the significance of the di- 
photon excess becomes substantially lower and agreement between theory and experiment 
can be reached at the ~ la level. This is achieved at the expense of slightly increasing the 
discrepancy of the H — >■ ZZ and H — )■ WW signals which, however, are affected by much 
larger experimental (mainly statistical) uncertainties. 

Let us start by discussing the two main theoretical uncertainties that enter into play in 
'^{99 ~^ H) and which have been discussed in detail by the LHC Higgs cross section working 
group (LHCHWG) [1] and in Refs. [SllIMO]: the scale and PDF+a, uncertainties. 

The perturbative QCD corrections to the gg ^ H cross section are so large, leading 
to a i^-factor of about 2.5, that it raises worries about the rate of convergence of the 
perturbative series. The possibility of still large higher order contributions beyond NNLO 
hence cannot be totally excluded. The effects of the unknown contributions are usually 
estimated from the variation of the cross section with the renormalisation /ir and factori- 
sation scales at which the process is evaluated. In the gg ^ H process, the median 
scale is taken to be /i_R = /iF = A^o = [S1II9] in order to absorb some of the soft-gluon 
resummation corrections. Indeed, a{gg — )■ H) calculated at NNLO with /iq = \Mh is then 
approximately the same (up to a few percent) as the resummed cross section at next-next- 
to-leading- logarithm (NNLL) with /iq = Mh [20]. However, the scale variation of the two 

^However, because Rzz seems to be in agreement with the SM and there is no sign of a new particle in 
direct searches at the LHC, the iJ — ?► 77 excess wih be particularly difficult to accommodate in relatively 
simple and/or well motivated SM extensions. One would probably have to resort to slightly "baroque" 
new constructions or scenarios to explain the excess. 

■^It might be noted that history could repeat itself: in the first Z — > (.^^~ events observed by both the 
UAl and UA2 Collaborations and which led to the discovery of the Z boson in 1983, a significant fraction 
were accompanied by additional photons jl5) . This triggered a plethora of papers proposing composite 
models of quarks, leptons and weak bosons, before the excess of photons died away with more statistics. 

^In particular, a significant fraction of the 77 events seems to come with two additional jets, while the 
predicted rate in the SM from vector boson fusion qq — > qqH and gg — > Hgg is expected to be small. 

^An example out of many for such a possibility is the pp — > bb production cross section at the Tevatron 
that had been first determined to be a factor of two to three larger than the QCD prediction, before higher 
order effects and various uncertainties were included; see Ref. |16| for a discussion. 

^An addendum with the complete analysis for Higgs production at y/s ~ 8 TeV has been added to the 
version of Ref. [5] submitted to the archives (,arXiv:1012.0530v 5). 
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cross sections is different. 

To estimate the scale uncertainty, the current convention is to vary the renormahsation 
and factorisation scales within the range iiq/k < ^r,^f < ^/^O; with the ratio of scales 
restricted to the range 1/k < ^p/f^R < f^- The choice k = 2 is usually adopted. For a 
Higgs mass Mh = 126 GeV, this leads to a scale uncertainty of Ao"^ ~^io% ~ '^^^ 

and Aa^ ^-9^5% = 8 TeV when the constraint |</iF//^_R<2is imposec£|. Slightly 

larger uncertainties occur if the scale is varied in a wider domain. If for instance, if one 
chooses K = 3, the scale variation would lead to a few percent more uncertainty. 

A second issue is related to the not yet entirely satisfactory determination of the parton 
distribution functions (PDFs) and in particular, the gluon densities. In addition to this, 
since a^'^i^gg ^ H) oc a1 and receives large contributions at 0{> a^), a modest change 
of as (which is possible as the average value [22] can be rather different from the ones 
obtained from deep-inelastic scattering data that are used in some of the PDFs) can lead 
to sizable change in the cross section value. There is a statistical method to estimate the 
PDF uncertainties by allowing a la (or more) excursion of the experimental data that 
is used to perform the global fits. In addition, the MSTW Collaboration [23] provides a 
scheme that allows for a combined evaluation of the PDF and uncertainties. 

To take into account this additional uncertainty and the spread in the predictions using 
the various NNLO PDF sets [23p^ . the PDF4LHC working group recommends [25] to take 
as a global PDF uncertainty the MSTW PDF+A^^Pos uncertainty at the 68% confidence 
level (CL) and multiply it by a factor of two. This procedure gives nearly the same answer 
as the one proposed in Refs. [TTJUB] in which one evaluates the combined 90% CL MSTW 
PDF+A^'^Ptts-l- A^^ctg uncertainty, where A^^as is for the error generated by the theoretical 
uncertainty on as estimated to be ~ 0.002 at NNLO [23]. For a{gg — H) at NNLO, one 
then finds a total PDF uncertainty of Aa^°^ ~ 9% at both ^/s = 7 and 8 TeV for a 126 
GeV Higgs boson. 

This discussion is illustrated in Fig. [1] where the gg ^ H inclusive cross section for 
a Mh = 126 GeV Higgs boson, evaluated at NNLO-QCD and including the electroweak 
corrections, is displayed as a function of the reduced scale yu//io with fi = fiR = and 
Ho = \Mh at ^/s = 8 TeV. The situation at ^/s = 7 TeV is very similar. One can see 
that indeed, for the usual choice fi/ fiQ = ^ (2) of scale variation, the cross section increases 
(decreases) by ~ 10%. If one is conservative and enlarges the domain of scale variation, 
one notices that (T{gg H) decreases monotonically with increasing /i//io but it has a 
plateau for about /x//xo ~ | where the cross section is maximal. 

In the left-hand side of Fig. [H the spread in the prediction for (r{gg — ?■ H) is shown 
for the six PDF sets that are available at NNLO [231121], including the pure 90%CL PDF 
uncertainty bands. One sees that the spread of the cross sections is rather significant, but 
most sets predict a rate that is smaller than the MSTW prediction except for NNPDF 
and HERAPDF. The right-hand side of the figure shows the PDFs uncertainty that the 
PDF4LHC group [25] recommends to retain, i.e. the MSTW PDF+A^^Pa, at 68%CL 
combined uncertainty multiplied by a factor of two, which is about ±9% in the entire /i/ /io 

^The uncertainty is a few % smaller if the two scales are equated, jiji ~ ^p, and one would obtain at 
NNLO Acr^ «^9'5% at = 8 TeV in accord with Ref. The small difference for the central value 

of the cross section, obtained in our case by using the latest version of the program HIGLU [21], is due to 
some refinements in the treatment of the electroweak corrections performed in Ref. [19] . 
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range. 

A critical issue is the way the scale and PDF uncertainties should be combined. As 
advocated by the LHCHWG [1], one should be conservative and add the two uncertainties 
linearly; this is equivalent of assuming that the PDF uncertainty is a pure theoretical un- 
certainty with a flat prioij^. For Mh ~ 126 GeV, this procedure leads to a total uncertainty 
of about /\scaie+PDF^ ^ ±20%. This is shown in the right-hand side of Fig. [1] 




Figure 1: a{gg H) at NNLO as a function of ^/^o for Mh = 126 GeV and ^ = 8 TeV. Left: 
the spread in the cross sections when the six NNLO PDF sets are used, normalised to the central 
MSTW cross section with jiQ = \Mh- Right: the relative PDE+a^ uncertainties in the MSTW 
case, when compared to the central value, as advocated by the LHCHWG group as well as the 
total uncertainty when the EFT uncertainty is linearly added. 

There is, however, a third source source of uncertainty which has not been accounted for 
by the LHCHWG [1] but has been discussed in Ref. [S1[I71[TB|. As the gluon-gluon fusion 
process, already at LO, occurs at the one-loop level with the additional complication of 
having to account for the finite mass of the loop particle, the NLO calculation is extremely 
complicated and the NNLO calculation a formidable task. Luckily, one can work in an 
effective field theory (EFT) approach in which the heavy quark in the loop is integrated 
out, making the calculation of the contributions beyond NLO possible. While this approach 
is justified for the dominant top quark contribution for Mh ^ 2mf [26], it is not valid for 
the 6-quark loop (and for the interference between the h- and the t-loops) and for those 
involving the electroweak gauge bosons [11] . The uncertainties induced by the use of the 
EFT approach at NNLO are estimated, from the NLO case in which both the exact and 
EFT calculations are available, to be of C(9%) for Mh ~ 126 GeV [5]. 

This uncertainty is of pure theoretical origin as it is due to an approximation in the 
calculation and has nothing to do with the scale uncertainty. Hence it should be added 
linearly to the scale+PDF uncertainty. For Mh ~ 126 GeV, this leads to a total uncertainty 
of A*°*(T ±30% on the NNLO gg H cross section when the scale is varied in the 
commonly adopted range | < /i//io < 2 as also shown in the right-hand side of Fig. [1] 

^Despite the fact that the Hessian method provides an error that is of probabilistic nature, it docs 
not account for the theoretical assumptions that enter into the PDF parametrisation. This theoretical 
uncertainty is reflected in the larger spread in the central values of the PDF predictions. 
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In the case of interest, i.e. for the normahsed cross section times branching ratios Rxx 
given by the ATLAS and CMS Collaborations, the story is not yet over and there are in 
fact two additional sources which lead to uncertainties or normalisation problems, albeit 
smaller than the ones discussed above: 

i) There are first uncertainties in the Higgs branching ratios. Indeed, while Higgs decays 
into leptons and gauge bosons are well under control, as mainly small electroweak effects 
are involved, the partial widths into quark pairs and gluons are subject to uncertainties. 
These are mainly due to the errors on input values of the bottom and charm quark masses 
and as, which then migrate to the other decays branching fractions [SlETj. For Mh ~ 126 
GeV, one would obtain ABR(77, WW, ZZ) ^ ±2% [5] when slightly smaller errors on the 
input masses and nic [27] compared to the PDG input values [22] are adopted. 

a) There is a slight problem with the overall normalisation of cr^gg — )■ H). The nor- 
malisation adopted by the experiments (and which comes from the LHCHWG) is the one 
obtained at NNLL [20] and not at NNLO [5l[T9]. Besides the fact that it is theoretically 
not entirely consistent to use the resummed result (as the PDFs are defined at NNLO 
and not NNLL), the resummation is available only for the inclusive rate and not for the 
cross sections when experimental cuts are incorporated and that are actually used by the 
experiments. It turns out that for Mh ~ 126 GeV, ct'^^'^ll -g ^ ^o/^ smaller than ct^^^^lo 
has a smaller scale dependence [20]. Hence, the gg H cross section might have been 
underestimated from the very beginning, albeit by only a small amount. 

These might affect the total rate and the uncertainties, increasing them by a few per- 
cent. If one also takes a non-dogmatic approach and increase the scale variation beyond 
the commonly chosen range 1 < fi/mo < 2, one might end up with a total theoretical 
uncertainty that is closer to A'^a ^ 40% than 30% for Mh ~ 126 GeV. 

We are now in a position to discuss the impact of this total uncertainty on the rates for 
Higgs production times decay branching ratios in the channels if — ?■ 77 and H ZZ that 
have been measured by the ATLAS and CMS Collaborations. A very important fact to note 
from the very beginning is that, in the experimental combination of different uncertainties, 
the theoretical uncertainty is not treated as a bias, as should be the cas^but, instead, as 
a nuisance parameter. Therefore, the scale and PDF uncertainties are not added linearly 
to the experimental uncertainties in contradiction with the LHCHWG recommendation, 
but are combined quadratically with the experimental statistical and systematical errors. 
As the latter are much larger than the scale and PDF uncertainties, at least 30% for the 
experimental errors and only 10% for the scale and 10% for the PDF uncertainties, the 
magic of statistics and the combination in quadrature makes that 

A*°V = V(A^-Pa)2 + (A'^(t)2 + (APDF^)2 ^ A"xp^ ^cxp^ ^ ^^^^ ^PDF^ 

This means that for A^^^a ^ 30%, which is the minimal experimental error, one would 

^Let us illustrate this important point by calculating a{gg — > H) in a consistent way but different from 
the one which gives the central value of Ref. |4] and which has been adopted as a normalisation by the 
ATLAS and CMS Collaborations. We choose to use the NNPDF2.1 set and evaluate the cross section at a 
scale fiR = fiF = fJ-o = jMh, which is within the range adopted for the scale uncertainty. We then obtain 
a NNLO cross section of a{gg H) ^ 22.9 pb for Mh = 126 GeV and 8 TeV. This value is w 20% 

larger than the reference value of <j{gg H) = 19.2 pb used by ATLAS and CMS. It is therefore clear 
that treating this theoretical uncertainty, which leads to the 20% change in the normalisation, as a mere 
nuisance can affect the conclusions in a very significant way. 
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obtain only A*°*(t ^ 33%. Hence, one can consider that, in practice, the above mentioned 
theoretical uncertainties are simply not reflected in the errors of eq. ([1]) for the ATLAS 
and CMS measurements of the rates in the different channels. 

For the comparison of theoretical predictions and the experimental measurements, it 
is more convenient to adopt the procedure advocated, for instance, in Ref. [TB], that is 
to ignore the theoretical uncertainties in the likelihood fit performed in the experimental 
analyses and simply confront the pure experimental error with the theoretical prediction 
that includes the theory uncertainty band. For the cross section in a given channel^ 
and because the experimental values obtained from the multivariate analyses have been 
cross-checked by a cut based analysis, it should be possible to use such a procedure. 

This is what is done in Fig. [2] for a Higgs mass of 126 GeV. First we combine the theory 
predictions for the gg^H cross section and their uncertainties at ^/s = 7 and 8 TeV and, 
because the integrated luminosity in both the 2011 and 2012 data samples is approximately 
the same, we simply perform an average. We then add all the theoretical uncertainties 
in two possible scenarios: i) only the scale and PDF uncertainties as advocated by the 
LHCHWG and which leads to Af^^^^y^^a ^ ±19%, ii) add linearly to the previous result 
the EFT uncertainty leading to a total of AJ;*Ypdf+eft'^~28%. 

In Fig. [21 the green and yellow bands represent these two possibilities for the total 
uncertainty. These bands are compared with the experimental measurements (which again 
we identify as a first approximation with the optimal value of the strength modifier fx) for 
the normalised production rates in the two channel^, R^^ and Rzz of the ATLAS and 
CMS Collaborations, as well as the ATLAS and CMS combination, given in eq. ([1]). 

It is clear that including the theoretical uncertainty helps to reduce the discrepancy 
between the experimental and theoretical values in the if — > 77 channel, while keeping 
the accord between the data and the SM prediction in the H ZZ channel. In the 
approach where the scale, PDF and EFT uncertainty are added linearly, one would obtain 
in the if — )■ 77 channel deviations with a significance of 0.7cr, 1.24cr and 1.3(T for the CMS, 
ATLAS and ATLAS©CMS resuhs, respectively. 

One should finally comment on the optimal value of the strength modifier when all 
channels are combined, /itot; given by the ATLAS and CMS Collaborations when the 2011 
and 2012 data are added: /itot = 1-2 ± 0.3 and /itot = 0.8 ± 0.22, respectively!^. If 
this parameter is to be viewed as a cross section measurement, it would mean that the 
gg ^ H is already "measured" to better than ^ 25% (since all channels analysed by the 
two collaborations are initiated by gg H, except for the H ^ bb channel for which 
the sensitivity is still rather low). This total error, which should be in principle largely 

-'^^ Although the situation here might be shghtly complicated as, in practice, one has to add to the cross 
section of the gg — > H process those of the vector boson fusion and Higgs-strahlung processes. However, 
because the gg — > H cross section is an order of magnitude larger than the summed cross section of the 
vector boson and Higgs-strahlung processes, the latter can be omitted in a first approximation. 

^^We do not include the gg ^ H ^ WW channel as first, it has not been yet analysed fully by the 
ATLAS Collaboration and second, the cross section in this case is broken into 0, 1 and 2 jet bins and this 
introduces an additional uncertainty due to the jet veto which can be significant [28j . 

^^As can be seen, there is an apparent deficit in the CMS Higgs cross sections despite the apparent 
excess in iJ — 77 rate. If the latter should turn out to be due to an upward statistical fluctuation and 
the global deficit should remain, the situation can be very easily accomodated by the 20-30% theoretical 
uncertainty in the normalisation of (j{gg — > H), that we have been discussing here. 
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Figure 2: The value of Rxx for the 77 and ZZ final states given by the ATLAS and CMS 
Collaborations, as well as their combination, compared to the theoretical uncertainty bands. 



due to the presently limited statistics, is of the same order of the theory uncertainty in 
the best case. We believe that this "paradox" will be resolved if the approach that we 
advocate, that is comparing the data for the cross sections including only the experimental 
uncertainties to the theoretical prediction with the uncertainty bands. 

In conclusion, we have first recalled that there are substantial theoretical uncertainties 
in the cross section for the dominant Higgs production channel at the LHC, gluon-gluon 
fusion, stemming from the scale dependence, the parton distribution functions and the 
use of an effective field theory approach to evaluate some higher order corrections. They 
are about 10% each and if they are combined according to the LHCHWG, they reach the 
level of 30% when the EFT uncertainty is also included. However, in the experimental 
analyses, these theoretical uncertainties in a{gg — )■ H) are treated as nuisance parameters 
rather than a bias. As they are still individually smaller than the experimental (statistical) 
errors, the net result is as if they had not been included in the total errors given by the 
ATLAS and CMS Collaborations. If the experimental results for the production cross 
sections times decay branching ratios in the various analysed channels are confronted with 
the theoretical prediction, including the theoretical uncertainty band, added linearly on top 
of the experimental error the discrepancy between the measurements and the prediction 
becomes smaller. This is particularly the case for cr[gg — )• H) x BR(iJ — t- 77), where the 
!^ 2cr discrepancy with the SM prediction reduces to the level of ^ la if the 30% theory 
uncertainty is properly considered. 
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